
* OPEN STATA OUTPUT FILE LOG *

log using "C:\Users\jp18390\Dropbox\DISCRIMINATION PROJECT (U.S. FEDERAL AGENCIES)\Age Discrimination Project\PAR R&R\Statistics\Age_Discrimination.APPENDIX 1.MISC SENSITIVITY CHECKS.06-02-2022.smcl" 
   
   

use "C:\Users\jp18390\Dropbox\DISCRIMINATION PROJECT (U.S. FEDERAL AGENCIES)\Age Discrimination Project\PAR R&R\Statistics\Age_Discrimination.Post-Estimation.06-02-2022.dta", replace 




**** 2022 PAR Data Replication [06/02/2022]: "STATUS-GROUP POWER DIFFERENTIALS AND AGE DISCRIMINATION WITHIN THE U.S. FEDERAL WORKFORCE" [KRAUSE & PARK] ****

 
 
**** APPENDIX SECTION 1: MISCELLANEOUS SENSITIVITY CHECKS *****
 
 
 
 
***  APPENDIX SECTION 1A: MODEL SPECIFICATIONS CONTROLLING FOR ONE-YEAR LEAD OF ORGANIZATIONAL JUSTICE COVARIATE [CONCERN: SIMULATANEITY BIAS OF THIS CONTROL COVARIATE WITH AGE DISCRIMINATION FORMAL COMPLAINTS MAY INADVERTENTLY BIAS THE MAIN COVARIATES OF INTEREST] **** 


***  APPENDIX SECTION 1B: MODEL SPECIFICATIONS DROPPING VETERANS AFFAIRS; DROPPING EXTREME CASES OF THE MAIN COVARIATES & DEPENDENT VARIABLE **** 

   
   
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
   


xtset a_id year, yearly

*
*   
*


** Create One-Year Ahead 'Lead' of Organizational Justice Covariate **

gen f1_orgjustice_sem = f1.orgjustice_sem




*** TABLE A1.1: APPENDIX SECTION 1A REGRESSION MODEL TABLE PREDICTING VARIATIONS IN AGE DISCRIMINATION FORMAL COMPLAINTS IN U.S. FEDERAL AGENCIES [REVERSE CAUSALITY TESTS: INCLUSION OF ONE-YEAR AHEAD "ORGANIZATIONAL JUSTICE" CONTROL COVARIATE] ***



** MODEL A1.1A: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio_40over_suplb  ratio_40over_nonsuplb   orgjustice_sem f1_orgjustice_sem   nonprof40over_tr_lb   politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10  ratio_55over_totallb, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio_40over_suplb  ratio_40over_nonsuplb    orgjustice_sem  f1_orgjustice_sem  nonprof40over_tr_lb   politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic
*
*

** COMPUTE IRR [EXPONENTIATED] MARGINAL EFFECT FOR PROCEDURAL FAIRNESS BASED ON BETWEEN INTERDECILE CHANGE IN COVARIATE **
lincom _b[B__orgjustice_sem]*0.2659021 - _b[B__orgjustice_sem]*-0.1089841, eform
*
*
*
*
*

**  MODEL A1.1B: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio40suplb_nonsuplb    orgjustice_sem  f1_orgjustice_sem  nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10  ratio_55over_totallb, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio40suplb_nonsuplb   orgjustice_sem  f1_orgjustice_sem   nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup  ratio_minsup_nonmsup  ratio_55over_totallb) 
*
estat ic
*
*
*


** COMPUTE IRR [EXPONENTIATED] MARGINAL EFFECT FOR PROCEDURAL FAIRNESS BASED ON BETWEEN INTERDECILE CHANGE IN COVARIATE **

lincom _b[B__orgjustice_sem]*0.2659021 - _b[B__orgjustice_sem]*-0.1089841, eform


*** THE ORGANIZATIONAL JUSTICE/PROCEDURAL FAIRNESS BETWEEN-EFFECT IS ATTENUATED IN MAGNITUDE AND ESTIMATED WITH IMPRECISION [LOSE SIGNIFICANCE] DUE TO COLLINEARITY WITH PRE-DETERMINED "LEAD" VERSION OF COVARIATE ***
*** NOTE: DATA IS NOT SUFFICIENTLY RICH ENOUGH TO DISCRIMINATE THESE CONTEMPORANEOUS AND LEAD ORGANIZATIONAL JUSTICE EFFECTS ON AGE DISCRIMINATION FORMAL COMPLAINTS *** 

correlate orgjustice_sem  f1_orgjustice_sem



*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************




*** TABLE A1.2: APPENDIX SECTION 1B: MODEL SPECIFICATIONS DROPPING VETERANS AFFAIRS; DROPPING EXTREME CASES OF THE MAIN COVARIATES & DEPENDENT VARIABLE ****  



** MODEL A1.2A: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio_40over_suplb  ratio_40over_nonsuplb   orgjustice_sem   nonprof40over_tr_lb   politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb  if a_id!=99 & a_id!=100 & a_id!=101, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio_40over_suplb  ratio_40over_nonsuplb    orgjustice_sem   nonprof40over_tr_lb   politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic
*
*
*
*
*
*
*



**  MODEL A1.2B: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***

xthybrid age_discrimination  ratio40suplb_nonsuplb    orgjustice_sem   nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10  ratio_55over_totallb if a_id!=99 & a_id!=100 & a_id!=101 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio40suplb_nonsuplb   orgjustice_sem    nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic
   
   

*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************

 


** MODEL A1.2C: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS [EXCLUDE VALUES OF DEPENDENT VARIABLE "AGE-DSICRIMINATION" EXCEEDING 99TH PERCENTILE VALUE [X > 386] ***

xthybrid age_discrimination  ratio_40over_suplb  ratio_40over_nonsuplb   orgjustice_sem   nonprof40over_tr_lb   politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb  if age_discrimination <= 386, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio_40over_suplb  ratio_40over_nonsuplb    orgjustice_sem   nonprof40over_tr_lb   politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic
*
*
*
*
*
*
*



**  MODEL A1.2D: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS [EXCLUDE VALUES OF DEPENDENT VARIABLE "AGE-DSICRIMINATION" EXCEEDING 99TH PERCENTILE VALUE [X > 386] ***

xthybrid age_discrimination  ratio40suplb_nonsuplb    orgjustice_sem   nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb  if age_discrimination <= 386, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio40suplb_nonsuplb   orgjustice_sem    nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic
   

*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************




** MODEL A1.2E: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS [EXCLUDE VALUES OF KEY COVARIATE "ratio_40over_suplb" EXCEEDING THE 1ST AND 99TH PERCENTILE VALUE] ***

xthybrid age_discrimination ratio_40over_suplb    ratio_40over_nonsuplb  orgjustice_sem    nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10  ratio_55over_totallb if ratio_40over_suplb>= 0.7095544 & ratio_40over_suplb<= 0.969697, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) ///
full use(ratio_40over_suplb    ratio_40over_nonsuplb   orgjustice_sem  nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb)

*
estat ic
*
*
*
*
*
*
*



**  MODEL A1.2F: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS [EXCLUDE VALUES OF KEY COVARIATE "ratio_40over_nonsuplb" EXCEEDING THE 1ST AND 99TH PERCENTILE VALUE]  ***

xthybrid age_discrimination ratio_40over_suplb  ratio_40over_nonsuplb   orgjustice_sem   nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if ratio_40over_nonsuplb>= 0.4495327 & ratio_40over_nonsuplb<= 0.8635558 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio_40over_suplb  ratio_40over_nonsuplb   orgjustice_sem    nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup  ratio_minsup_nonmsup  ratio_55over_totallb) 
*
estat ic
   

*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************
*****************************************************************************************************************************************************************************************




** MODEL A1.2G: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS [EXCLUDE VALUES OF KEY COVARIATE "ratio40suplb_nonsuplb" EXCEEDING THE 1ST AND 99TH PERCENTILE VALUE] ***

xthybrid age_discrimination  ratio40suplb_nonsuplb    orgjustice_sem   nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if ratio40suplb_nonsuplb  >= 1.028541 &  ratio40suplb_nonsuplb  <= 1.834515, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
use(ratio40suplb_nonsuplb   orgjustice_sem    nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 
*
estat ic
*
*
*
*
*
*
*

   

log close
